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Abstract. We propose a novel method for sorting events of multiparticle production according to the 
azimuthal anisotropy of their momentum distribution. Although the method is quite general, we advocate 
its use in analysis of ultra-relativistic heavy-ion collisions where large number of hadrons is produced. The 
advantage of our method is that it can automatically sort out samples of events with histograms that 
indicate similar distributions of hadrons. It takes into account the whole measured histograms with all 
orders of anisotropy instead of a specific observable (e.g. V 2 , 52 )- R can be used for more exclusive 

experimental studies of flow anisotropies which are then more easily compared to theoretical calculations. 
It may also be useful in the construction of mixed-events background for correlation studies as it allows to 
select events with similar momentum distribution. 

PACS. 25.75.-q Relativistic heavy-ion collisions - 25.75.Gz Particle correlations and fluctuations - 
02.50.Ng Distribution theory and Monte Carlo studies 


1 Introduction 

Hot matter which is created in ultrarelativistic heavy- 
ion collisions expands very fast in both longitudinal and 
transverse directions [T1I2I3I1]. The expansion is always 
anisotropic. This is true even in most central collisions, 
where one would expect symmetry in azimuthal angle due 
to circular shape of the initial overlap region [5]. In non¬ 
central collisions, the overlap of the two colliding nuclei 
has elliptic shape and thus naturally second-order (and 
higher even orders) anisotropy in the fireball expansion 
builds up. On top of this—at any centrality—energy de¬ 
position in the interactions of the incoming partons fluc¬ 
tuates and this leads to all orders of anisotropy in the 
transverse expansion of the fireball. Therefore, even within 
carefully selected centrality class flow anisotropies vary 
from event to event mi. 

Even more degrees of freedom are offered by colli¬ 
sions of deformed nuclei, like uranium. There, the initial 
anisotropy will also depend on the way how the colliding 
nuclei are oriented. 

This makes the comparison of experimental data to 
theoretical simulations more complicated, because one has 
to take into account that every event starts with different 
initial conditions and evolves differently. Simulations are 
compared to data in order to pin down the properties of 
the matter which is being modelled. Initial conditions are 
unknown, however, although recent hydrodynamic results 
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indicate that their fluctuations can be directly mapped 
onto measured fluctuations of hadron distributions m- 
Comparison of theory to data must be done with great 
care so that the spectra of theoretical and experimental 
fluctuations match each other. Experimentally, events are 
distributed into centrality classes according to multiplic¬ 
ity. There is a problem with this procedure on the side 
of theory if very narrow centrality class is demanded, e.g. 
ultra-central collisions corresponding to 0-0.2% centrality. 
Events with the same multiplicity may evolve from ini¬ 
tial states with different impact parameter^ Moreover, 
all those events would differ by the quantum fluctuations 
in initial energy and momentum deposition [10] . There¬ 
fore, events from a class selected according to multiplicity- 
based centrality may have evolved from different initial 
conditions and experienced quite different evolution his- 


^ This can be seen, e.g., in Fig. 2 of [9], where the procedure 
is explained that is used by ALICE collaboration to deter¬ 
mine centrality. Within the used Monte Carlo Glauber model 
if the impact parameter is fixed, then the number of partici¬ 
pants may still fluctuate. On the other hand, from the overlap 
of centrality classes in Vpart histogram it is clear that fixed 
Vpart corresponds to an interval of impact parameters. As the 
experimental multiplicity is determined from multiplicity in a 
chosen detector, there is yet another source of fluctuation that 
comes from the uncertainty between multiplicity and Vpart- In 
summary, if we fix the multiplicity—even in perfectly spheri¬ 
cal nuclei like Pb or An—there is still some freedom for the 
impact parameter to fluctuate which we estimate of the order 
1 fm. Even more fluctuations can be expected in collisions of 
non-spherical nuclei, like U. 
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tory. It would be useful if there was a more selective method 
to choose events that are more likely to have evolved sim¬ 
ilarly. In addition to this, in collisions of deformed nu¬ 
clei the multiplicity itself is definitely not sufficient selec¬ 
tion criterion since the same multiplicity may result from 
events with very different initial orientations of the col¬ 
liding nuclei and thus very different flow patterns. Again, 
the situation calls for a selection method like the one pre¬ 
sented here. 

It would thus be advantageous if one could select a 
collection of events among all measured events (which may 
or may not belong to the same centrality class) which 
show very similar distribution of the produced hadrons. 
For such events one can assume that they also evolved 
similarly. 

The word “similar” when talking about distributions 
or histograms can be understood in layman’s terms so that 
they have similar shapes when you rotate them appropri¬ 
ately. There is a possibility how to quantify this with the 
help of the distance measure in the Kolmogorov-Smirnov 
test, but this will not be used here because it does not au¬ 
tomatically provide a way for sorting. We rather use the 
Bayesian framework where we basically ask the question: 
how would the events be grouped, based on the shapes of 
their histograms. 

A method aiming for such event selection has been 
proposed in m and is commonly referred to as “Event 
Shape Engineering”. It usually employs the size of the 
flow vector (/„ defined on a selection of hadrons from a 
given event (usually referred to as subevent) as 

Ms Ms \ 

^cos(n<).,), ^ sin(n(;i,) j , (1) 

i=i i=i ) 

where Ms is the multiplicity of the subevent and (pi are 
the azimuthal angles of the individual hadrons from the 
subevent. Events are then selected based on their values 
of qn (in most cases q 2 ). However, only the subevents not 
used in q^ determination can be used in further studies in 
order to avoid bias. 

It is not clear, moreover, whether selecting events ac¬ 
cording to the value of provides the best possible selec¬ 
tion method aimed at collecting similar events. It actually 
may not be the case, as we demonstrate below. 

In this paper we propose the use of a novel method for 
comparing, sorting and selecting events according to sim¬ 
ilarity with each other. The method is adopted from [H 
|T3] and was previously used in a different context [T3] . Its 
uniqueness consists in not fixing a single observable which 
would then be used for sorting of the events. It rather 
compares complete histograms (e.g. azimuthal angle dis¬ 
tributions) of individual events and it sorts the events in 
such a way that events with similar histogram shapes end 
up close together. After such a sorting has been performed 
one simply selects similar events just by choosing a group 
of events which follow each other in the created series. 

On the selected groups of events one could measure 
various observables (f„’s, g„’s, radial flow, temperature, 

...) which should fluctuate much less than in the whole 
measured event sample. 


A natural application of the method is in construction 
of correlation functions. There, one often needs a reference 
distribution which is constructed via the mixed events 
technique. If events used in mixing are different, this may 
introduce unwanted artificial effects into the correlation 
function. Therefore, the mixed events sample is always 
constructed with aligned event planes. (For application at 
intermediate energy nuclear collisions see e.g. [mns]). 

With the help of the proposed method it would be in¬ 
teresting to perform femtoscopic studies where oscillations 
of radii in azimuthal angle in both second and third order 
together at the same time should be visible. 

We comment more on the applications in the Outlook 
section. 

In the next Section we shall explain the method and in 
Section[3]we illustrate its use on Monte Carlo data from a 
toy model. The method is applied on more realistic Monte 
Carlo data generated by a transport model in Section 
We conclude in Section and give an Outlook about pos¬ 
sible applications of the method and its next development. 


2 The method 

Suppose that we have a sample consisting of a large num¬ 
ber of events. Initially, we can sort and divide those events 
into N percentiles according to the value of a single ob¬ 
servable which can be measured in every event. This can 
be the value of < 72 , V 2 , multiplicity or any other observable. 
Generally, we shall refer to this observable as Q. For the 
sake of clarity let us explain the method with a particular 
choice: choose = 10 event bins (deciles) and the observ¬ 
able according to which we sort data is Q = 92 = | 92 |- Let 
us stress at this point that the method is universal and 
will not depend on the choice of N and Q. 

In each of the event bins we can now produce the his¬ 
togram of hadron distribution in azimuthal angle summed 
over all events in the event bin. There is no physics in how 
the two nuclei are oriented when they collide, and so we 
have a free choice of how to rotate individual events before 
adding them to the angle histogram. For the introduction 
of the method let us align each event according to the sec¬ 
ond order event plane. Note, however, that the choice of 
alignment is a sensitive issue which we shall discuss later. 

Thus each event is characterized by the bin record of 
its distribution in azimuthal angl^ {ui} and belongs to 
one of the N event bins which we number with (We 
shall use Latin letters for angle bins and Greek letters for 
event bins.) 

2.1 Basic relations 

Imagine that we take a random event from our big sample 
and ask an unbiased observer, in which event bin he or she 

^ We denote Ui the number of particles in angle bin i and 
the whole record of an event is referred to with the help of 
braces. Thus summation of all angle bin entries gives the event 
multiplicityy),. Ui = M. 





Renata Kopecna, Boris Tomasik: Event shape sorting 


3 


thinks that this event belongs. More specifically, we can 
ask the question in the framework of Bayesian probability: 
What is the probabilitjj^(p,|{ni}) that an event with bin 
record {rii} belongs to the event bin /r? 

This probability will be evaluated with the help of 
Bayes’ theorem 


When the formula ([^ is inserted into ([^ and Q we 
obtain the practically usable relation from "vraich the large 
factorials drop out 






( 8 ) 


( 2 ) 

where P{A\B) is the conditional probability of the event 
A given event B. The probability of event B can be de¬ 
termined 

P{B) = Y,PiB\A)P{A), (3) 

A 

where the sum runs over all possible events A. Definitions 
of symbols P{B\A) and P{A) are analogical. 

With the help of Bayes’ theorem we can express 


With the help of this conditional probability we can 
determine for an event with angle bin record {rii} its mean 
event hin number 

Jl = ^nP{n\{ni}). (9) 

u 

2.2 The algorithm 

Now we are able to describe the algorithm which is used 
for sorting of the events. 


P{fi\{ni}) 


P{{n^}\fi)P{fj:) 

PiM) 


( 4 ) 


Here, P({ni}|/i) is the probability that if one randomly 
draws an event from the distribution function given by 
average histogram of event bin /i, the result will be the 
bin record {ui}. The prior P{fi) = 1/N takes the value of 
0.1 now. The denominator contains the overall probability 
of drawing the event {ui} from any of the event bins. It 
can be determined according to eq. ([^ 

N 

P{{ni}) = P{{n,}\iy)P{iy) . (5) 

1^=1 


1. First, all events are sorted according to the observable 
Q. 

2. Events are divided into N event bins according to cur¬ 
rent sorting. 

3. For each event and all event bins the probability 
P{^\{ni}) is determined that the event with record 
{ui} belongs to event bin /r. The mean event number 
p, is calculated for each event according to 

4. Events are sorted again according to their values of /x. 

5. Events are divided into N event bins according to cur¬ 
rent sorting. 

6. If the new sorting changed assignment of any events 
into event bins, the algorithm returns to step 3. Oth¬ 
erwise it converged. 


The advantage of using the latter formula comes from the 
fact that we are able to determine P({ni}|z/) for each bin 
record {ni\ and every event bin v 

p(K)i.)=Min Slip- (6) 

t 

Here M is event multiplicity, the product goes over all 
angle bins i and P{i\v) is the conditional probability that 
random particle falls into angle bin i given that the event 
to which it belongs stems from event bin v. It can be 
determined when we take the number of particles from all 
events in event bin v falling into angle bin i. This number 
is divided by the total number of all particles from events 
in event bin n, denoted M^, 


® Obviously, the event must belong to one of the event bins 
so that the probabilities must be normalised 

N 

^P(p|{ni}) = 1. 


The construction of the algorithm is such that once con¬ 
verged, the events are sorted so that those ending up close 
to each other are characterized by similar angular his¬ 
tograms. This is the best possible experimental approach 
to the selection of events that have undergone similar evo¬ 
lution. This is as good working definition of what similar 
events are as it can be. 

In the present formulation the average histograms are 
more strongly determined by high-multiplicity events. Note, 
however, that the algorithm is independent of event mul¬ 
tiplicity. Of course, if in particular physics analysis certain 
multiplicity is demanded, one can easily use it for event 
selection and then use the present method for more refined 
selection of the events. 

The final sorting of events also does not depend on 
the initial sorting. Hence, this method can also be used 
for a judgment if the particular observable Q, e.g. (72 or 
U 2 , is a good measure for selecting similar events. If the 
initial ordering according to Q is correlated with the final 
ordering, then the observable Q is good for this purpose. 
This may not always be the case, as we will show later. 

Although the final result does not depend on initial 
ordering, a good initial ordering can lead to faster conver¬ 
gence of the algorithm. 

There is a caveat, however, which has been discussed 
only shortly so far. The algorithm works well for sorting 
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Table 1. Parameters used for generating multiplicity depen¬ 
dent Vn- 


n 

an X 10^ 

bn X lO"* 

Cn 

1 

0 

0.01667 

-0.000680 

2 

-7.099 

20.06 

0.07874 

3 

-2.083 

6.658 

0.04236 

4 

-96.38 

2.621 

0.04897 

5 

-71.76 

2.236 

0.01673 


histograms in observables that are not periodic, e.g. rapid¬ 
ity. However, the azimuthal angle of particle momentum 
is a periodic observable, i.e. we can always arbitrarily ro¬ 
tate the events. Practically, two angular histograms may 
be almost identical when they are both aligned properly, 
but might appear quite different for the proposed algo¬ 
rithm if they are rotated in random directions. Thus the 
way how events are aligned initially plays a crucial role 
and we observed that it strongly biases the final sorting 
of the events. At the moment we do not have recommen¬ 
dation for an automatic algorithm which would align the 
events in the best way. Instead, we tested a few reasonable 
choices for the initial alignment in our toy model studies. 
With the simple toy model, quite naturally that initial ro¬ 
tation aligning the second-order event planes 1^2 will lead 
to sorting characterized by V 2 - Analogically, with align¬ 
ing the third-order event planes ips, final sorting is given 
by V 3 . In the next Section we present results which take 
into account a combination of both these orderings. That 
feature was, however, weaker in AMPT-generated events. 


3 Illustration of the method 

3.1 Elliptic flow 

Let us first demonstrate the action of the sorting algorithm 
on a simple case of events with only first- and second- 
order anisotropic flow. We generated the azimuthal angles 
of pions from the distribution 

-^ 2 ( 0 ) = (1 + 2ui cos(^ - V'l) + 2u2 cos(2(0 - ■02))) ■ 

Ztt 

.. 

The parameters vi and V 2 depend on the multiplicity of 
the event M as 

Vn = anM"^ + b„M + Cn , (H) 



0 0.2 0.4 0.6 0.8 1 

Initial sorting 



V2 

Fig. 1. Top: correlation of the resulting sorting variable p with 
the initial ordering. Every dot represents one event. Bottom: 
correlation of p. with V 2 determined for each event via the event 
plane method. 

of second-order event planes defined from the generated 
Monte Carlo data for each event via 

M 

^ ( 12 ) 


and parameters a„, bn, and c„ for each n can be found in 
Table [l] They have been de term ined from experim enta l 
data as we discuss in Section 3.2 In addition to eq. (II), 
flow anisotropy parameters Vn are Gaussian-smeared with 
a width of 0.25. For each event, the directions of event 
planes ipi and rp 2 are random and uncorrelated. 

We generated 5000 events with multiplicities between 
300 and 3000. Directed flow is practically negligible and 
the dominant anisotropy is second-order. Thus the most 
reasonable choice of initial event rotation is the alignment 


To show the power of the method we first ordered the 
events fully randomly. The algorithm converged after 65 
iterations. In Fig. top one can see that the initial or¬ 
dering indeed has nothing to do with the final ordering of 
events. In the bottom panel of that Figure we show the 
correlation of final sorting variable p, with the value of V 2 
determined in each event via the event-plane method (re¬ 
sults from cummulant method are practically identical). 
One can clearly see that the ordering is given by the el¬ 
liptic anisotropy of the particle distribution. 
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Note that when we started the sorting algorithm with 
initial ordering according to the value of q 2 in each event, 
the final correlation was just a mirror image of Fig. [^bot¬ 
tom. Small values of V 2 corresponded to small ft and high 
values of V 2 to high values of ft. This illustrates that the al¬ 
gorithm always converges to a sorting of events according 
to their similarity but the direction how they are ordered 
along p, may be different. 

There is a congestion of events seen at p = 9 and 10 
in the upper panel of Fig. and a step in the same place 
in the lower panel of that figure. This actually shows that 
the sorting was very clear for these two event bins. The 
congestion is made out of events for which when evaluating 
p according to eq. (§, the sum has one clearly dominant 
term. The probabilities that the event might belong to 
other event bins are very small. 

For illustration, in Fig. [^we show average histograms 
of the events in the individual event bins. We see that the 
relative amplitude of the second order variation (i.e. the 
amplitude divided by the mean value of the bins) decreases 
from event bin 1 to event bin 10. 

The method is thus able to distinguish different shapes 
of hadron distributions and sort events according to them. 
In this simple case, this might not look like a big advantage 
as we could have sorted the events simply by measuring V 2 ■ 
Therefore, we proceed with a more complicated example 
where the algorithm demonstrates its full power. 


3.2 Anisotropic flow 


We parametrised the dependence of w„’s on multiplicity 
through a fit to data from ALICE and ATLAS collabo¬ 
rations [TTIIIH]. Coefficients Vn for n = 1 through 5 are 
parametrised according to eq. 0 with parameters sum¬ 
marised in Table [T] 

Then we generated 5000 events with multiplicities be¬ 
tween 300 and 3000 pions and angular distributions ac¬ 
cording to 


^*5 (</>)= ^ ( 1 + ^ 2VnCOs{n{4>- lijn)) 


(13) 


The n„’s for every event are set by eq. (11) and then 
smeared with Gaussian distribution with the width of 
0.25. The phases ipn are selected from uniform distribution 
and are not correlated with each other. 

Now we have to address the question how to rotate the 
events so that the comparison of individual events to the 
event bin histograms yields the most reasonable sorting. 
There are two symmetries at play here: rotational symme¬ 
try and parity symmetry. We can rotate an event around 
the collision axis, and we can also flip it so that we get its 
mirror image. We have observed that both these symme¬ 
tries influence the result. 

The two dominant components of flow anisotropy are 
second and third order. Hence, in our tests we focused on 
the corresponding event planes. If events are all aligned 
into the direction of second-order event plane, the algo¬ 
rithm becomes sensitive to the second-order anisotropy 


and to large extent ignores the third order. Analogically, 
alignment according to third-order event plane enhances 
the sensitivity to third order anisotropy. Thus the result¬ 
ing sorting is rather sensitive to this choice. Furthermore, 
in both cases the algorithm distinguishes events which 
look like mirror images of each other (i.e. have opposite 
parity). This must be taken into account when designing 
the sorting algorithm. 

The shape of hadron distribution is never solely deter¬ 
mined by V 2 or by V 3 . It rather follows from their com¬ 
bination which is different in every event due to varying 
phase difference of second and third order event plane. 
In order to take into account both these components of 
anisotropy we rotated all events according to the angle bi¬ 
sector between ip 2 and '03- We denote its azimuthal angle 
02-3- Also, in order to take care of the parity symmetry 
the events were oriented so that 02 is less than 7r/2 away 
from 02-3 counterclockwise. 

The initial sorting of the events was random and the 
algorithm converged after 121 iterations. 

From Fig. |^it is clearly seen that in this case sorting 
of the events is neither determined by q 2 , nor by W 2 , nor 
by V 3 . Higher order terms also do not play a big role at all. 
The event shape is complex and results from an interplay 
of all its simple characteristics. The message of the Figure 
is that 92 may not be a good variable to select events ac¬ 
cording to their shape because panel (a) shows that it is 
not correlated with the overall shape of the event as soon 
as more flow harmonics are involved. Note that in our toy 
model there is neither correlation between the flow har¬ 
monics nor between the event planes of different order. 
This may not be so in real events and then the correla¬ 
tion between sorting variable p and some of the measured 
quantities may appear. What we show is thus rather an 
extreme case. It calls, however, for attention: the overall 
shape of an event and thus the evolution dynamics running 
in that event cannot be simplified into a single measured 
variable. There might be a counterargument that variables 
like e.g. 92 which are used in Event Shape Engineering are 
proved to be good at event sorting, because events with 
different values of 92 show different values of other mea¬ 
sured quantities. However, in addition to this, our method 
optimizes the sorting so that events which are placed close 
to each other share as many event shape characteristics as 
possible. Eigure shows that such a sorting may not be 
connected with any of the commonly used variables. 

In Fig. we again observe that the resulting values of 
p are grouped around integers. This actually means that 
the assignment of the events into event bins is very clear, 
because in determining the value of p from eq. the 
probability P{p\{ni\) is (close to) one for certain pL and 
(nearly) zero elsewhere. The event then clearly belongs to 
the event bin /r. 

When we tried to start the sorting algorithm with ini¬ 
tial ordering according to the value of 92 it failed to con¬ 
verge within reasonable time (5000 iterations). This hap¬ 
pens sometimes when unfavourable initial ordering is used. 
On the other hand, in other occasions we have checked 
that also with different initial ordering of events the algo- 
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4i [rad] 

Fig. 2. Average histograms of the azimuthal angles for event bins 1-10, with event bins indicated in the panels. Events with 
only 112 and Vi. 


rithm converged to identical final sorting. The only differ¬ 
ence could arise from the feature that the algorithm does 
not follow any specific condition in which direction the 
events should be sorted. Thus different initial orderings 
may end up in mutually reversed final orderings. 

In Fig. we show the resulting average angular his¬ 
tograms after the sorting. The events indeed differ by their 
shape, not just only by the value of one of the flow har¬ 
monics. Since the second and third-order anisotropies are 
dominant and due to the initial rotation of the events, 
higher order harmonics are washed out by averaging over 
the event bins and not seen in the Figure. 

A question may appear to what extent the algorithm 
would sort the events even if they would be drawn from 
the same distribution. It would still try to place closer 
together events with similar histograms, and further apart 
those events with more different histograms. Then—if in 
doubt—one can test the hypothesis that events are drawn 
from the same distribution e.g. with the method proposed 
in m or with the help of Kolmogorov-Smirnov test. 

4 Application to AMPT events 

After we have established and tuned the sorting algorithm, 
we now use it in a more realistic setting with events gen¬ 
erated by the AMPT model ppj . 

The model was used as commonly distributed with two 
modifications which are recommended for realistic simu¬ 


lation of Pb-|-Pb collisions at the LHC collision energy 
■^/snn = 2.76 TeV: the parton screening mass was re-set 
to 2.097 fm“^ and the string melting was turned on in or¬ 
der to avoid the underestimation of partonic effects [5D]. 

We have generated 2000 events which correspond to 
the 0-20% centrality class. On the generated particles we 
have applied rapidity cut \y\ < 1 in order to roughly sim¬ 
ulate the acceptance of central tracking detectors. For the 
first rough analysis we have taken all charged hadrons and 
run the sorting algorithm on them. 

The event shapes are dominated by the second order 
anisotropy, but it is not the only feature that determines 
the shape of the fireball. When running our sorting al¬ 
gorithm, we have tried three different initial alignments: 
according to ip 2 , '03 and 02 - 3 - In all three cases only V 2 
and q 2 of the events show any pattern of correlation with 
the final sorting. This is shown in Fig.[^ Surprisingly, even 
in the case of 03 initial alignment we have found no cor¬ 
relation of final ordering with V 3 . On the other hand, as 
seen in Fig.[^ there is a pattern that shows that also here 
the sorting of events is strongly influenced by the second- 
order anisotropy. It will be interesting in the future to see 
if this dominance of second-order anisotropy survives also 
in other centrality classes, particularly in more exclusively 
selected central event. 

It is also interesting to see that in realistic simulation, 
alignment with respect to 02 does not automatically lead 
to such a clear correlation of V 2 and p, as it was the case 
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Fig. 3. Correlation of various observables with final sorting variable p. Simulated are events with anisotropies up to 5th order 
and initial rotation is according to ip 2 - 3 - Correlation with a) V 2 , b) va, c) < 72 , d) event multiplicity. 


with our toy model. This is seen in the upper two panels of 
Fig. [5] Such a correlation exists for the two or three event 
bins with highest /r’s. There, higher /r corresponds to a 
higher value of V 2 - However, for /2 below 7 there seems 
to be no correlation between the ordering of the event 
and V 2 - There is slightly more correlation in case weighted 
V 2 evaluation, as seen in Fig. We would like to un¬ 
derstand where the difference between the events shapes 
comes from, which forces them into different event bins. 
To this aim, we show in Fig. the average histograms 
in individual event bins after the sorting procedure has 
converged. In event bins 8 , 9, and 10 the gradual growth 
of V 2 is evident. In the other event bins the histograms 
show a more complicated structure, where higher order 
anisotropies also give an important contribution. We re¬ 
call, however, that no clear correlation of any higher Vn 
with the obtained fi is observed. 


The next suspected cause of the difference of events is 
the relative angle between the directions of 92 and 93 . We 
thus studied the correlation between p, and the relative 
angle. In most cases, no correlation was observed. There 
is a hint of correlation, though, in case that the events are 
aligned according to ^^ 3 . This is shown in Fig. Indeed, 
there seems to be a slight correspondence between the as¬ 
signment to an event bin and the angle between q 2 and 
93 . The interpretation is at hand, that for the event shape 
the relative phase is decisive. Unfortunately, this correla¬ 
tion is gone when we align the events differently. Thus, 
there are hints that the relative angle is important, but 
the assignment of an event to an event bin appears to be 
given by more complicated interplay of various individual 
features. 
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ij) [rad] 

Fig. 4. Average histograms of the azimuthal angles for event bins 1-10, with event bins indicated in the panels. Events with 
anisotropies up to 5th order. 


5 Conclusions and outlook 

It is very useful to have a method able to sort events 
in such a way that it is possible to select those with 
very similar momentum distribution. One can then assume 
that they must have undergone similar evolution and this 
makes it possible to study the dynamics of hot expanding 
matter more exclusively. 

We tested this method on artificial events generated 
with AMPT. It showed that dividing the events into classes 
according one selected variable, usually q 2 , does not really 
correspond to selecting event with the same shapes. It still 
seems that the main role in determining the event shape— 
at least for the studied centrality class 0-20%—is played 
by the second-order anisotropy. Nevertheless, other fea¬ 
tures are important as well. In one case we could identify 
the difference ip 2 — ' 4’3 to co-determine the assignment to 
event bins, but this observation is not universal for any 
initial event alignment and any way of evaluation of q^s. 

In addition to the explanation of the method of Event 
Shape Sorting, we thus gave a first superficial study of 
realistic events with the proposed method. The study of 
event shapes generated by AMPT for various centralities 
which would include thorough analysis of all features that 
influence the final shape and their physics interpretation 
would go beyond the scope of this work and we plan to 
publish it in a separate paper. 

As was mentioned already in the Introduction, the pre¬ 
sented method allows to select groups of events with simi¬ 


lar momentum distributions. If used in data analysis, one 
can then measure various quantities on such events and 
study how they are related to the event shape. 

The method might even allow to get as close to single¬ 
event femtoscopy as possible. First, when doing femtoscopy 
with a single event one would run into difficulties with the 
uncorrelated reference distribution which is usually con¬ 
structed through the event mixing technique. Event shape 
sorting could provide a selection of events with similar mo¬ 
mentum distributions which would make a suitable sample 
for event mixing. Second, the statistics in a single event 
would be too small to perform a 3D analysis. However, 
one could take a sample of events with similar momentum 
distributions and reasonably expect that they also have 
the same sizes and undergo the same dynamics. Then, 
one could analyse the correlation function integrated over 
the whole selected event sample. The feasibility of such 
studies will be investigated in the future. 

An interesting application appears to be the classifi¬ 
cation of events from U-|-U collisions. Due to deforma¬ 
tion of the colliding nuclei one expects large fluctuating 
anisotropies of the transverse flow. Surprisingly, a pre¬ 
liminary study of azimuthally sensitive correlation radii 
showed no dependence on the value of q 2 [H]. The latter 
was employed for the selection of events with different final 
state anisotropy. A short inspection of our Fig. would 
suggest that this is no surprise at all! As soon as there are 
other harmonic components of the anisotropy, the shape 
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Fig. 5. Correlation of the final sorting variable ^ with V 2 of individual events for events generated by AMPT. a) Initial alignment 
of the events according to ' 02 , b) initial alignment according to 02 and V 2 evaluated withpt weight, c) initial alignment according 
to 03 , d) initial alignment according to 11 : 2 - 3 ■ 


of the events is more complex. The “proper” partition of 
events into various classes by the type of anisotropy should 
be done differently. Our algorithm can do such a proper 
classification. 

Let us also comment again on the interesting though 
perhaps academic question, how the proposed algorithm 
would proceed if all events were indeed generated from 
the same underlying probability distribution and the only 
differences between them would be due to statistical fluc¬ 
tuations. The algorithm would be sensitive to the differ¬ 
ences whatever their cause might be. Thus it would sort 
the events so that typical fluctuations within one event 
class would be below the normal statistical ones. Such a 
situation could be detected with the help of standard sta¬ 
tistical tests, like e.g. the Kolmogorov-Smirnov test. In 
real data we do not expect this to happen, however. 


Also, there are still technical issues which require some 
discussion and will be addressed in the future. Most im¬ 
portant is the ambiguity if initial rotation of the events 
for which we do not yet have optimised rules. Another is 
the rather high requirement on CPU time for even mod¬ 
erately large event samples. Note however, that we have 
not tried any fancy computational optimisation of the al¬ 
gorithm so far, and hence we would expect some room for 
improvement here. 

In fact, we also work on a well optimized routine that 
can readily be taken and applied directly in data analysis. 
Integration into standard packages like ROOT or HistFit- 
ter [22] will be addressed, as well. 

In spite of this, we believe that the Event Shape Sort¬ 
ing is worthwhile to apply in real data analysis and carries 
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Fig. 6 . Average histograms of the azimuthal angles for event bins 1-10, from events generated with AMPT. Initial alignment 
of events according to ip 2 - 
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